Spin wave contribution to the nuclear spin-lattice relaxation in triplet superconductors 
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We discuss collective spin wave excitations in triplet superconductors with an easy axis anisotropy 
for the order parameter. Using a microscopic model for interacting electrons we estimate the fre- 
quency of such excitations in Bechgaard salts and ruthenate superconductors to be one and twenty 
GHz respectively. We introduce an effective bosonic model to describe spin-wave excitations and 
calculate their contribution to the nuclear spin lattice relaxation rate. We find that in the exper- 
imentally relevant regime of temperatures, this mechanism leads to the power law scaling of 1/Tj 
with temperature. For two and three dimensional systems the scaling exponents are three and five 
respectively. We discuss experimental manifestations of the spin wave mechanism of the nuclear 
spin lattice relaxation. 



Nuclear magnetic resonance (NMR) is a powerful tool 
for analyzing ordered electron states in solids. NMR 
analysis has been successfully applied to study magnetic 
insulators (see e.g. and references therein) as well 
as several classes of unconventional superconductors, in- 
cluding high Tc cuprates Q , heavy fermion materials Q , 
ruthenates 0, and organic superconductors 0. In par- 
ticular, NMR experiments have been useful for analyzing 
the symmetry of the SC order parameter |6( and for clar- 
ifying the structure of the phase diagram in systems with 
competing orders Q. 

A common feature of the NMR experiments in certain 
families of triplet superconductors (TSC) is the power 
law temperature dependence of the nuclear spin lattice 
relaxation rate. Bechgaard salts H, ruthenates and 
heavy fermion materials @ showed 1/Ti - T 3 at low 
temperatures and for small magnetic fields. Such behav- 
ior has usually been interpreted as a signature of nodes in 
the quasiparticle gap on the Fermi surface. Indeed, point 
and line nodes should lead to T 5 and T 3 scaling of 1/Ti 
respectively. In several cases, however, we have reasons 
to doubt the presence of nodes in the TSC order param- 
eter. For example, in quasi one dimensional Bechgaard 
salts the natural order parameter has different signs on 
the two sheets of the Fermi surface and no nodal points 
|lfj|. Also in ruthenates, the order parameter that is con- 
sistent with spontaneous time reversal breaking |llj | and 
the quasi two dimensional nature of these materials cor- 
responds to a constant quasiparticle gap on the entire 
Fermi surface [l^. I n this paper we consider a mecha- 
nism of the nuclear spin lattice relaxation that is not due 
to Bogoliubov quasiparticles but due to collective spin 
wave (SW) excitations of the TSC order parameter. We 
demonstrate that in the experimentally relevant regime 
of temperatures, this mechanism also leads to the power 
law scaling of 1/Ti with temperature. 

Our starting point is the Moriya relation 0] for the 
nuclear spin lattice relaxation rate 



J_ = 2g^\Al 
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Here A describes the strength of hyperfine interactions 
between nuclear spins and conduction electrons, g^ is 
a gyromagnetic ratio of the nucleus, g e s is an effective 
gyromagnetic ratio of conducting electrons, ps is a Bohr 
magneton, X±jj(<7> u iv) is the imaginary part of the trans- 
verse (i.e. perpendicular to the magnetic field) electron 
spin susceptibility taken at the nuclear Larmor frequency 
ojn- In the case of a perfect spin SU(2) symmetry, lin- 
early dispersing SW excitations exist down to arbitrar- 
ily small energies. In real materials there is always spin 
anisotropy which gives rise to a finite gap for spin ex- 
citations, ujq. Below, we estimate the value of ojq to be 
tens of millikclvin for Bechgaard salts and hundreds of 
millikclvin for the ruthenates. This is much larger than 
the nuclear Larmor frequency, lon, but smaller than the 
typical temperature used in experiments. When ujq is 
much larger than ujn, creation and annihilation of in- 
dividual SWs does not affect x"(wjv). However, there 
is a contribution due to the scattering of thermally ex- 
cited SW excitations. Let p(E) be the density of states 
for SW excitations and n(E) = (exp(E/k B T) - l)" 1 be 
the Bosc distribution function. From the second order 
perturbation theory we have x'L^n) ~ / p{E)p(E + 
u>ti)[n(E) — n(E + uj]y)]dE. The characteristic energy 
scale in this integral is set by the temperature T. Since 
T > wo we can assume linear dispersion of SW exci- 
tations and take p(E) ~ £/ x , where d is the num- 
ber of spatial dimensions. Using u>n <C T, we have 
x" z (un) ~ lo n J E 2d - 2 (-dn/dE)dE ~ T 2d ~ 2 . Com- 
bining this result with the Moriya relation (JIJ we obtain 
1/Ti ~ T 2d_1 . This simple analysis does not take into 
account coherence factors in the expression for \" ■ Be- 
low we demonstrate that coherence factors do not modify 
the scaling exponent of the nuclear spin lattice relaxation 
rate in TSC. This is in contrast to antiferromagnetic sys- 
tems that also have linearly spin waves, but in which 
coherence factors contribute an additional 1/T 2 factor 
to 1/Ti Q. 

For a detailed analysis of the nuclear spin lattice relax- 
ation rate we introduce an effective model that captures 
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the essence of triplet superconductivity and allows us to 
analyze collective excitations. A simple picture of the 
TSC state corresponds to binding electrons into Cooper 
pairs with spin one and momentum zero and Bose con- 
densing such pairs. In an effective bosonic model one can 
neglect details of the orbital nature of Cooper pairs and 
consider them as "elementary" particles. Interactions are 
important for the correct description of collective excita- 
tions, thus we are led to the Hubbard type model for spin 
one bosons [Til ] 



n = -t 



(13)° 
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Here a ia creates a boson on site i with spin a = 
{ — 1,0,1}. Operators rii and Si describe the number 
of atoms and the total spin on site i: Uj = a i<7 ai a , 

= So-.cr' a \<7 T <J°' a i<y' i wner e T aa > are the usual spin 
matrices for spin one particles. The first term in the 
Hamiltonian describes tunneling of Cooper pairs be- 
tween neighboring lattice sites i and j. An important 
aspect of the model is the presence of two types of 
interaction terms. The third term in J3J) depends only on 
the number of particles on each site and is the same as for 
spinless bosons. The fourth term in (0) gives spin depen- 
dence to the interaction (without breaking the spin SU(2) 
symmetry) and is a novel feature of spinful Cooper pairs. 
The sign of U2 determines the difference between unitary 
(U2 > 0) and nonunitary (U2 < 0) triplet superconduc- 
tors. It is generally believed that triplet pairing between 
fermions leads to unitary Cooper pairs |l2(. Thus, from 
now on we assume that U2 is positive. For concreteness, 
we consider a d dimensional (d = 2,3) hypercubic lat- 
tice. Our results, however, do not depend on the precise 
lattice structure. The second term in equation |J2j in- 
troduces easy axis anisotropy for the order parameter by 
making the condensation of ao to be energetically favor- 
able. A state with finite (ao) corresponds to the unitary 
state of Cooper pairs with the d vector pointing along 
the z axis (see eq. (|HJ for the definition of the d vector). 

In the mean field approximation, we take (ao) = \&o 
with |* | 2 = (p + zt + 6r)/Uo and z — 2d being the 
coordination number. Without loss of generality, we 
can take \I/o to be real. Fluctuations in the phase of 
ao correspond to the density (Bogoliubov) mode. To 
find SW excitations we need to consider a± operators. 
In the Hamiltonian (J2J we replace ao by its expecta- 
tion value, take the terms quadratic in a±, and perform 
the Bogoliubov rotation a ia = ^= Ha^^ a k+ = 

v kl% + + u kl^_j:_, a -k- = Uk ~?l + + v kl-k-- Here N is 
the total number of lattice sites, u^ + v 2 = (£& + A)/Ek, 
2u k v k = -A/E k , & = -2tY^ d n=1 cos{k n b) + U^l-^. In 
these equations b is the lattice constant and A = t^^n- 



We obtain the diagonalized spin- wave Hamiltonian H.q = 
E fc E k (j k+lk+ + 7 |_ 7fc _) with El = e k + 2A£ fc . Op- 

erators jI ± create SW excitations with S z = ±1. In 
the long wavelength limit we find E 2 = ojq + v 2 k 2 with 



= 5r 2 + 2ASr and 



2t U 2 ^ 2 a 2 



We need to calculate the electron spin susceptibility in 
the direction perpendicular to the applied magnetic field. 
Let 6 be the angle between the z axis (i.e. the direction 
of the d vector) and the direction of the magnetic field 
(see Fig. Pi . We have 



X±h = sin 2 6 Xzz + (1 + cos 2 9) Xx 



(3) 



It is easy to see that Xxx and Xzz are given by the corre- 
lation function of SWs 

Xxx(q,u) = 2{g cff ^ B ^ ) 2 (u q + v q ) 2 D %{l,u) 

a/3 

Xzz(q,u) = (5cffMs) 2 J U af}(k,k + q) x 



D^ a (k + q,n + Lj)DUk,Q) 



(4) 



where we introduced the Nambu-Gorkov type nota- 
tions D*< A (q,uj) = J dt e^*0(±t)(T ^ qa {t)^{0)), with 

^ qa = {l q +,l- q -} T , and U a f}{k,k') = 8 a p{v k v k , - 
u k u k i) + (1 - 6 a f/)(u k v k , - v k u k i). 

Note that there is a qualitative difference in calculat- 
ing Xzz and Xxx- A non-uniform magnetic field in the z 
direction can scatter the existing thermally excited SWs. 
Thus we find finite imaginary part of Xzz at small fre- 
quencies by considering the quadratic Bogoliubov Hamil- 
tonian. We obtain 

x '^N,q) = 2(g c ^ B ) 2 ^b d J^L 



x n(E k )(n(E k+q ) + l)S(E k - E k+q ). 



(5) 



By contrast, a non-uniform magnetic field in the x direc- 
tion can only create or annihilate SW excitations. How- 
ever, energies of these excitations cannot be smaller than 
luq. Hence, if we limit ourselves to the quadratic Hamil- 
tonian for SWs, we find that Xxx IS identically zero for 
frequencies smaller than ujq at any temperature. To get 
finite Xxx at small frequencies we need to consider inter- 
actions between SWs. Taking quartic terms in equation 
(J2J and using definitions of SW operators, we obtain the 
interaction terms between SW excitations. These can 
be used to calculate self-energies for SW excitations as 
shown in Fig [2 We find [3 

x'Uq,^) = W B *o) 2 ^(t/o + 3L/ 2 ) 2 fe 2!i 



d d h f d d k 2 E 2 kl + E- 



E k3 



(2ir) d J (2ir) d E kl E k2 E k3 
x (n(E kl ) + l){n(E k2 ) + l)n(E k3 ) 
x 6(E kl +E k2 - E k3 ). 



(6) 
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where k% — k\ + &2- We emphasize that equations (0 
and JfjJ) apply only in the low frequency limit wjv <C ojo 
which is relevant for experiments. 

Expressions (J5J and JBJ show that in equation © con- 
tributions to 1/Ti due to x'L and x"a; scale as T 2d ~ x and 
rp3d-2 respectively. For two and three dimensional sys- 
tems the x'zz contribution dominates (one can check that 
this conclusion remains when we include prefactors). At 
first glance this result appears surprising. Firstly, the 
real static susceptibility is finite in the direction perpen- 
dicular to the d vector but is zero along it. Secondly, in 
the case of full SO (3) symmetry, creation and annihila- 
tion of individual SWs contribute to x'xx and Xy V at small 
frequencies, but have no effect on x'zz- A crucial part of 
our analysis is the existence of a spin gap loq which is 
much larger than the nuclear Larmor frequency ujn- In 
this case x' '(ww) does not have contributions due to cre- 
ation or annihilation of individual SWs. For x'zz we take 
thermally excited SWs and scatter them by the magnetic 
field. For x'xx we also need thermally excited SWs but 
in addition we must rely on interactions between them. 
SWs are pseudo-Goldstone modes. At low energies inter- 
actions between them are suppressed. This gives rise to 
the smallness of x'xx relative to x'L ■ 

To summarize, for two and three dimensional systems, 
we find that as long as is not anomalously small, the 
SW contribution to the nuclear spin lattice relaxation 
rate is given by 

Here S d = 2n d / 2 /T(d/2) is the surface area of a unit 
sphere, and B n are Bernoulli numbers, B 2 = 1/6, B4 = 
— 1/30. We remind the readers that equation applies 
when T ^> ujq. At low temperatures, T <C ojq, we expect 
1 /Ti to start decreasing exponentially, reflecting the ex- 
ponential suppression in the number of thermally excited 
SW excitations. Equation (JJJ) also predicts that the nu- 
clear spin lattice relaxation rate should be very sensitive 
to the direction of the magnetic field. We note, however, 
that this argument is valid only for magnetic fields that 
are smaller than the so-called spin-flop magnetic field, 
-ffflop- In magnetic fields larger than -fffl op , the order pa- 
rameter d will always be perpendicular to the direction 
of the applied field [15| . 

It is useful to compare contributions to 1/Ti from 
magnons to the one from quasiparticles. For concrete- 
ness we consider a quasi two dimensional system with 
the TSC order parameter d — Aozk x /kp, that has a line 
of nodes along the z axis. The quasiparticle contribu- 
tion to the NMR relaxation rate in such a state [l6j is 
given by 1/Ti qp = tt 2 \A\ 2 g^b 2 m 2 T 3 /6A 2 ,. For compar- 
ison, we take equation (0 and use the BCS expressions 
for the velocity of SW excitations (see below). We find 
Timag/Tiqp = 7t 2 Aq/£|,. In a typical superconductor, 



the value of the quasiparticle gap is much smaller than 
the electron Fermi energy. Therefore, when we have both 
gapless quasiparticles and SWs, the quasiparticle contri- 
bution will strongly dominate. Only when the quasipar- 
ticles are fully gapped out do the magnons provide the 
dominant contribution to 1/Ti. 

Now we outline the key steps of the analysis that al- 
lowed us to estimate the characteristic frequency of SW 
excitations in TSC. We consider a phenomenological BCS 
type model for interacting electrons 

« = E e *4«rCto-E V a _1 (?)4(?)da(9) 
kcr aq 

d{(q) = \^Va{q)f k cl +q/2a {ia 2 a a ) aP c^ k+q/2l3 {%) 

kafi 

Here d a (q) are Fourier components of the TSC order pa- 
rameter in the direction a (a = x,y,z), fk is an orbital 
wavefunction (e.g. fk — sign(fc) [lflj), and ct CT are elec- 
tron creation operators. Parameters V a describe elec- 
tron interactions in the p-wave channel. A homogeneous 
unitary triplet superconducting state is the ground state 
of the Hamiltonian ijSJ. Assuming easy axis anisotropy 
with V z > V x = V v , we find (d z (q = 0)} = A and 
(d x ) — (d y ) — 0. To find SW excitations we consider fluc- 
tuations of the TSC order parameter d x (r, t) = d* (r, t) . 
Integrating out fermions gives an effective action for 
SW excitations [3 Seff{4 } = J uq D x (q,uj) \ d x (q,u) | 2 , 
where q and lo are the wavevector and the real frequency 
of the SW. For the spherically symmetric Fermi surface in 
d dimensions, we have D x {u> 7 q) = jNo(lu 2 — ui 2 + v 2 q 2 ). 
Here is the density of states at the Fermi energy, 
v s = Vp/d, vf is the Fermi velocity, and 

c 2 = 4A 2 A (r 1 (^ 1 -K" 1 ) (9) 

Zeroes of the function D x correspond to SW excitations. 
In the case of spin symmetric interactions with V z =V X , 
SWs are Goldstone modes of broken spin symmetry and 
are gapless. They have linear dispersion with the velocity 
v s . In the case of the easy axis anisotropy, SWs have a 
gap oj . 

As a concrete example, we consider the TSC state in 
Bechgaard salts. We assume that the spin anisotropic 
part of the interaction in this phase is the same as in 
the antiferromagnetic state of this material A7i an ; s = 
SJz^^j} SfSj . Here (ij) corresponds to the nearest 
neighbor sites and the spin z axis points along the crys- 
tallographic b' axis. The antiferromagnetic resonance 
experiments 01 suggest SJ Z = 0.01K. We can ex- 
press A7i an i s in the form similar to equation JSJ) with 
-SV Z = 5V X = SV y = \5J z v . Here v is the volume 
of the unit cell. Assuming that the anisotropic term is 
a small correction to the spin symmetric interaction, the 
total interaction entering equation JSJl is V a — V + SV a . 
The value of V can be estimated from the BCS equation 



for the transition temperature T c — 1.14 WBose -1 ^ ^- 
Here wbos is the characteristic frequency of bosons pro- 
viding electron pairing. Combining all expressions, we 
find 



CJ 



= 2{SJ z N Q v ) 1/2 log( 



1.14w BO s 



(10) 



Taking SJ Z = 0.01K, N = 2 • 10 33 erg 
1000K, T c = 1.4 if, A = 2.5 K, and v a = 360 1 3 we 
obtain ujq around one GHz. There is a simple qualitative 
argument that supports our result that lvq for Bechgaard 
salts lies in the GHz range. This argument relies on a 
comparison of SW resonances in the antiferromagnetic 
and superconducting states of these materials. In quasi 
one-dimensional systems microscopic descriptions of the 
two states are similar and BCS type models and our anal- 
ysis of SW excitations can be used in the antiferromag- 
netic phase as well. From the discussion above, we expect 
that the ratio of SW energies in the two phases is approx- 
imately proportional to the ratio of AF and TSC tran- 
sition temperatures, i.e. around ten. In the AF phase 
the SW resonance frequency was measured in the tens of 
GHz range 0. Hence, in the superconducting state it 
should be a factor of ten smaller, which brings us into 
the GHz range. 

Similar analysis can be done for the TSC state in 
S^RuOi. Sigrist and coworkers (see e.g. Ref [isj ) 
showed that spin-orbit coupling in these materials leads 
to a difference in the transition temperature of various 
order parameters of the order of two percent. From the 
BCS expression for Tc we have ST C /T C = SV/N a V 2 . Tak- 
ing Ao = 2.6K, from equation jnj we estimate u>q to be 
around twenty GHz. 

Si^RuOi is a quasi two dimensional material. Bech- 
gaard salts have a mixed dimensionality [lil |20| with 
easy, intermediate, and hard directions in transport. 
Hence, they may also exhibit properties of a quasi two di- 
mensional system with respect to SW excitations. Both 
the ruthenates and Bechgaard salt superconductors ex- 
hibit the T 3 dependence of 1/Ti at low temperature. 
Therefore these materials should be good candidates for 
experimental investigation of the SW mechanism of nu- 
clear spin lattice relaxation discussed in this paper. 

In summary, we studied spin wave excitations in triplet 
superconductors with the easy axis spin anisotropy. We 
derived an explicit expression for the energy of such ex- 
citations and used it to estimate the spin wave energy in 
Bechgaard salts and ruthenate superconductors. We con- 
sidered the effects of spin wave excitations on the nuclear 
spin lattice relaxation rate. We showed that in the ex- 
perimentally relevant regime of temperatures, 1/Ti has a 
power law scaling with temperature, including the T 3 de- 
pendence for two dimensional systems. We showed that 
the spin wave mechanism predicts a dramatic decrease 
of 1/Ti for temperatures lower then the energy of spin 
wave excitations and leads to dependence of 1/Ti on the 



_1 cm 3 , wbos 




12 t + ^ y 



21 y 
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FIG. 2: Second order self-energy diagrams that gives rise to 
scattering of spin waves. For details see [l5|. 



direction of the applied magnetic field. 
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